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Semiclassical catastrophes in the dynamics of a quantum 
rotor (molecule) driven by a strong time- varying field are con- 
sidered. We show that for strong enough fields, a sharp peak 
in the rotor angular distribution can be achieved via time- 
domain focusing phenomenon, followed by the formation of 
angular rainbows and glory-like angular structures. Several 
scenarios leading to the enhanced angular squeezing are pro- 
posed that use specially designed and optimized sequences 
of pulses. The predicted effects can be observed in many 
processes, ranging from molecular alignment (orientation) by 
laser fields to heavy-ion collisions, and the squeezing of cold 
atoms in a pulsed optical lattice. 

PACS numbers: 42.50.-p, 42.50.Vk, 32.80.Pj, 33.80.-b 



I. INTRODUCTION 

Driven rotor is a standard model in classical and quan- 
tum nonlinear dynamics studies [Q. Increasing interest 
in the problem has arisen because of recent atom op- 
tics realization of the quantum (5-kickcd rotor |^,^, and 
novel experiments on molecule orientation (alignment) by 
strong laser fields [^-Q. A strong enough laser field cre- 
ates the so-called pendular molecular states § |l^ that 
are hybrids of field- free rotor eigenstates. By adiabati- 
cally turning on the laser field, it is possible to trap a 
molecule in the ground pendular state, thus leading to 
molecular alignment. The only way to reach a consider- 
able degree of alignment in this approach is by increas- 
ing the intensity of the field. However, many applica- 
tions may require only a transient molecular alignment 
(orientation), where the molecular angular distribution 
becomes extremely squeezed at some predetermined mo- 
ment of time. It is well known that a physically related 
problem of squeezed states generation in a harmonic sys- 
tem may be solved by a proper time-modulation of the 
driving force (parametric resonance excitation). Behav- 
ior of a rotor in general, strong, time-varying fields is 
a much less-studied problem, although it is understood 
that the long-persisting beats in the molecular angular 
distribution may be induced by short laser pulses [pd|-p^ . 

In the present paper, we analyze generic features in 
the dynamics of a quantum rotor driven by strong pulses, 
and present a strategy for cfhcient squeezing of the rotor 
angular distribution by a sequence of pulses of moder- 
ate intensity. The results of our research are related to 
a number of physical processes, ranging from molecular 



alignment (orientation) by laser fields to heavy-ion colli- 
sions, and the trapping of cold atoms by a standing light 
wave. The paper is organized as following. Section II 
discusses semiclassical catastrophes in the dynamics of 
a driven quantum rotor with the help of a simple two- 
dimensional model ||l^. Section III applies these generic 
results to a thermal system of cold atoms driven by a 
pulsed optical lattice Section IV focuses on 3D ef- 
fects important in the processes of molecular alignment 
and orientation by strong laser pulses pO|. 



II. TWO-DIMENSIONAL ROTOR 

For the sake of clarity, we start with the simplest 
model of a two-dimensional rotor described by the time- 
dependent Hamiltonian 



(1) 



where L is the angular momentum operator, and / is the 
momentum of inertia of the rotor. This model contains 
most of the physics we want to present here. Some impor- 
tant additional effects that appear in 3D will be discussed 
below in Section IV. For a "linear 2D molecule" having 
a permanent dipole moment fi, and driven by a linearly 
polarized field, the interaction potential is 



V{e,t) = -fi£{t) cos{9) 



(2) 



where £{t) is the field amplitude (i. e., of a half-cycle 
pulse), and 9 is the polar angle between the molecular 
axis and the field direction. In the absence of interac- 
tion with a permanent dipole moment, the external field 
may couple with induced anisotropic molecular polariza- 
tion. This interaction (being averaged over fast optical 
oscillations) may result in the interaction potential pro- 
portional to cos^(6l) (see Eq.(|l^) in Section IV) Al- 
though these two forms of V{9,t) may lead to different 
physical consequences (i. e., orientation vs alignment), 
the effects we will present are more or less insensitive to 
the choice of interaction. Therefore, we prefer to start 
with the more simple situation described by Eq.(^). 

By introducing dimensionless time r = th/I, and in- 
teraction strength e — iJ,£{t)I /h^ , the Hamiltonian can 
be written as 
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The wave function of the system can be expanded in the 
eigenfunctions of a free rotor 
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In the absence of the field, the wave function takes the 
form 



, +00 
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Despite a simple form of Eq.(^), the wave function ex- 
hibits extremely rich space-time dynamics. In particu- 
lar, it shows periodic behavior in time with the period 
Trev — 47r (full rcvival) and a number of fractional re- 
vivals at T = p/sTrev {p and s are mutually prime num- 
bers) . An analytical solution valid for a general time- 
dependent field is unknown even for this simplest model. 
Much effort has been devoted to the case of extremely 
short field pulses (^— kicks) (see, e. g. and references 
therein). In general, as a result of a single kick applied 
to the rotor at t = Tk, the coefhcients c„ transform as 



Cn{Tk +0) ^ 
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'Jn-m{P)Cm{Tk-0), (4) 



where 



P = 



e{T)dT, 



and Jn{P) is the Bessel function of the nth order. The re- 
sult of multiple kicks applied at different times can be ob- 
tained by combining transformations (^ after each kick 
with a free evolution according to Eq.(^ between the 
kicks. If the kicks are applied periodically to the system 
with the period T^ev, the system does not show chaotic 
behavior, and the energy accumulates quadratically with 
time (the so-called "quantum resonance" 24|). It is, 
therefore, quite natural to examine potential accumula- 
tion of angular squeezing of the rotor wave function under 
the "quantum resonance" excitation. In this case, be- 
cause of the exact quantum revivals at the free-evolution 
stages, the effect of N kicks of a magnitude P is equiva- 
lent to the action of a single strong pulse of strength NP 
(see, e. g. [^). In Figure 1, we show numerically calcu- 
lated time evolution of the probability density |5'(6',t)P 
after a relatively strong kick of a magnitude P = 85 
applied at r = 0. Initially the rotor was in the ground 
s-state (c„(0) = (5,io). For the chosen values of r, several 
distinct phenomena can be seen in these plots. First of 
all, the wave function shows an extreme narrowing in the 
region of small 9 after some delay following the kick [Fig. 
1 (b)] . The physics of this effect may be understood with 



the help of the following semiclassical arguments. Con- 
sider an ensemble of randomly oriented classical rotors 
subject to a kick. The angular velocity of a rotor located 
at the angle 9 is 



Lo{e) = -Psin[9) 



(5) 



just after the kick, assuming negligible initial velocity. 
For rotors from the region of small << 1, the acquired 
velocity is linearly proportional to their initial angle, so 
that all of them arrive at the focal point ~ at the 
same time 



(6) 



This phenomenon is quite similar to the focusing of light 
rays by a thin optical lens. For P >> 1, the shape of 
the distribution at the focusing time tj is dictated by 
the aberration mechanism (deviation of the cos{9) po- 
tential from the parabolic one), and it is P-independent. 
We consider the orientation factor O —< 1 — cos(6') > 
(where angular brackets mean averaging over the state 
of rotor) as a measure of the rotor orientation. For large 
enough P, the time-dependent orientation factor (for the 
initial s-state) may be easily estimated by averaging over 
the initially uniform classical ensemble of rotors having 
the velocity distribution of Eq. (|): 0(r) = 1 - Ji(Pr). 
Here Ji{x) is Bessel function of the first order. The min- 
imal value O « 0.418 of the orientation factor is, in fact, 
achieved in the post-focusing regime, at t « 1.84t/. As 
seen in Fig. 1 (c) , a new phenomenon can be observed in 
the angular probability distribution just after the focus- 
ing. Sharp singular-like features are formed in the distri- 
bution, which are moving with time. Each of these fea- 
tures has a typical asymmetric shape, with pronounced 
oscillations on one side and an abrupt drop down on the 
other side. Again, the origin of this effect can be traced 
in the time evolution of a classical ensemble of initially 
motionless rotors. 

After a kick applied at t = 0, the motion of the rotors 
is described by 



61 = 6*0 - Psin(6lo)r (mod 27r) 



(7) 



where 6*0 is the initial angle. For t < t/ Eq.(0) represents 
a one-to-one mapping 9{9q) [see Fig. 2 (a)]. At r = r/ 
the curve 9{9o) touches the horizontal axis [Fig. 2 (b)]. 
At T > Tf the angle becomes a multi- valued function of 
9 [Figs. 2 (c), (d)]. The classical time-dependent angular 
distribution function of the ensemble is given by 



\d9/d9^\ 



(8) 



The summation in Eq. (|^) is performed over all possible 
branches of the function 9o{9) defined by Eq.(^. It fol- 
lows immediately from Eq. (^) that even for a smooth ini- 
tial distribution, f{9, r) may exhibit a singular behavior 
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near the angles where d9/d9'^ 0. The quantum nature 
of the rotor motion replaces the classical singularities by 
sharp maxima in the probability distribution with the 
Airy-like shape typical to rainbow phenomena. Indeed, 
this effect is similar to the formation of caustics in the 
wave optics p5[ |, and rainbow- type scattering in optics 
and quantum mechanics p^-pst- We should stress, how- 
ever, that the long-time asymptotic regimes are radically 
different for the classical and truly quantum motion of 
the rotor. Thus, contrary to the classical limit, in which 
the caustics exist forever, they gradually disappear in the 
quantum case because of the overall decay of the initial 
rotational wave packet. On even longer time scale, an- 
other quantum phenomenon can be seen, namely revivals 
and fractional revivals of the initial classical-like motion. 
Figs. 1 (e)-(i) show several examples of fractional foci and 
rainbows in the angular distribution, which is a purely 
quantum effect. 

As we have demonstrated, a mere application of 
(5— kicks at the condition of "quantum resonance" does 
not lead to accumulated angular squeezing, and the orien- 
tation is saturated at some finite asymptotic level. Here 
we suggest an excitation scheme that exhibits the de- 
sired accumulation property. As previously mentioned, 
the wave function of the rotor reaches the state of the 
maximal orientation (i. e., minimal O value) after a cer- 
tain delay Ati following the application of the first kick 
at T = Ti =0. We suggest to apply the second kick at 
T2 = Ari. Immediately after the second kick, the sys- 
tem will keep the same probability density distribution. 
On the other hand, t = T2 will no longer be a station- 
ary point for 0{t) =< 1 — cos(6') > (r). The orientation 
factor 0{t), and its derivative are continuous and peri- 
odic functions of time in the course of a free evolution. 
Therefore, 0(t) will reach a new minimum at some point 
T2 + At2 in the interval [t2, T2 + Tr^v]- Clearly, the new 
minimal value of the orientation factor is smaller than 
the previous one. By continuing this way, we will apply 
short kicks at iterative time instants t^+i — + At^;. 
By construction of this pulse sequence, the squeezing ef- 
fect will accumulate with time, in contrast to the "quan- 
tum resonance" excitation. This is demonstrated by 
Fig. 3, which shows calculated sequences {Arfe} and 
{< 1 — cos(6') > (t^)} for a rotor initially in the s-state 
and being kicked by pulses with P = 3. The logarithm 
of the orientation factor gradually decreases, without any 
sign of saturation. 

At the stage of a well-developed squeezing (O << 1), 
the cos{d) -potential may be approximated by a parabolic 
one. It can be easily shown (both classically and quan- 
tum mechanically) that in this limit our strategy pro- 
vides the following recurrent relationships for the time 
intervals Ar^, successive values of the angular variance 
Uk —< >k~ '2'Ok and normalized variance of the an- 
gular velocity Wk —< {~id/dd)^ >k /-P^ : 



Uk 



Uk + Wk 



Uk+1 = Uk - 
Wk+1 = Wk 



Uk 
Uk 



Wk 



(9) 



For large fc, the last two finite-difference equations may 
be replaced by a system of coupled differential equations. 
The latter has an exact solution providing the following 
asymptotics: < 9^ >fecx: 1/Vk and At^ cx 1/k. This 
result is in good agreement with the numerically observed 
power-laws behavior of the graphs 3 (a) and 3 (b), and 
it describes correctly their slopes at k >> 1. We note, 
that in contrast to the wave optics (in which the size 
of the focal spot is diffraction limited), our system may 
be, in principle, " unlimitedly" squeezed in angle. We 
also note that a quasi-periodic sequence of kicks applied 
at Tk+i — Tk + Arfc -|- Trev provides the same squeezing 
scenario for a quantum rotor. The introduction of the 
Tret,-shift between pulses may be useful in the practical 
realizations of the scenario to avoid the overlap between 
short excitation pulses of a finite duration. 



III. SQUEEZING OF ATOMS IN A PULSED 
OPTICAL LATTICE 



In the present Section we apply the above results to 
another well-known system : cold atoms interacting with 
a pulsed optical lattice |19|. Optical lattices are periodic 
potentials for neutral atoms induced by standing light 
waves formed by counter-propagating laser beams. When 
these waves are detuned from any atomic resonance, the 
ac Stark shift of the ground atomic state leads to a con- 
servative periodic potential with spatial period A/2, half 
the laser wavelength (for a review, see, e. g. |^). Such 
lattices present a convenient model systems for solid state 
physics and nonlinear dynamics studies. In contrast to 
traditional solid state objects, the parameters of optical 
lattices (i. e. lattice constant, potential well depth, etc.) 
are easily controllable. Many fine phenomena that were 
long discussed in solid state physics, have been recently 
observed in corresponding atom optics systems. Bloch 
oscillations |Q, or the Wannier-Stark ladder are 
only few examples to mention. Time-modulation of the 
frequency and intensity of the constituent laser beams 
provide tools for effective modeling of numerous time- 
dependent nonlinear phenomena. Since the initial pro- 
posal and first pioneering experiments on atom op- 
tics realization of the (5-kicked quantum rotor cold 
atoms in optical lattices provide also new grounds for 
experiments on quantum chaos. 

We describe atoms as two-level systems with transition 
frequency u>o, interacting with a standing optical wave 
that is linearly polarized and has the frequency of w;. 
If the detuning A; = ujq — uji is large compared to the 
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relaxation rate of the excited atomic state, the internal 
structure of the atoms can be neglected and they can be 
regarded as point-like particles. In this approximation, 
the Hamiltonian for the atomic motion is 

2 

H{x,p,t)^ ^-V{t) COS {2kix), (10) 

where m is the atomic mass, fc; = wj/c is the wave num- 
ber of the standing wave. The depth of the potential 
produced by the standing wave is V{t) — fiil{t)^ /8Ai, 
where il{t) = 2\dzE{t)\/fi is the Rabi frequency, d is the 
atomic dipole moment and E(t) is the time-dependent 
amplitude of the light field. In the case of rather short 
laser pulses, this Hamiltonian corresponds to that of the 
two-dimensional 5-kicked rotor considered in the previous 
Section, with 9 = 2kix. In accordance with the previous 
discussion, many aspects of the dynamics of this system 
can be explained with semiclassical arguments. Here we 
provide results for a classical description of atoms in a 
pulsed optical lattice, with thermal effects taken into ac- 
count. 

We performed a Monte-Carlo simulation of the dynam- 
ics of an ensemble of (5-kicked particles which were ini- 
tially uniformly distributed in space and had a thermal 
momentum distribution. Figure 4 shows the spatial dis- 
tribution of atoms at different times after a single kick. 
In the upper row, the initial temperature of the ensem- 
ble is zero, and we can observe focusing [Fig. 4(a)] and 
formation of caustics [Fig. 4(b)] as it was described in 
Section II. In Fig. 4(c) and (d), the average thermal en- 
ergy of the initial ensemble is comparable with the energy 
supplied by a kick. As a result, instead of sharp peaks 
in the spatial distribution (like in Figure 4 (a) and (b)), 
we observe some broader spatial structures that are still 
reminiscent of the focusing phenomenon and the rainbow 
effect. 

We examined the accumulative squeezing approach in 
application to the above system. The angular localiza- 
tion factor O =< 1 — cos {2kix) > introduced in Section 
II describes now the spatial width of atomic groups lo- 
calized in the minima of the light-induced potential. In 
the classical limit, the localization factor after applying 
n kicks is given by 

0(t„) = 1 - — / dujo d0op{ujo)cos0n, (11) 

where 0„ is determined by 

0n=On-l+ATn(^UJo~'j2sme}j. (12) 

Here, luq — 2kiVo, where vq is the initial velocity of an 
atom, and picoo) describes the thermal distribution of 
the velocities. Figure 5 displays the amount of spatial 
squeezing for the series of 100 kicks for different initial 



temperatures. It can be seen, that after the first few kicks 
the localization factor demonstrates a negative power de- 
pendence as a function of the kick number (a straight line 
in the double logarithmic scale). Although, the system 
demonstrates a reduced squeezing for higher initial tem- 
peratures, the slope of all of the curves in Figure 5 is 
the same after the first several kicks, in full agreement 
with the general arguments of the previous Section (see 
Eqs.(^. Therefore, the accumulative squeezing scenario 
may be an effective and regular strategy for atomic lo- 
calization even at finite temperatures. 

However, this does not mean that the accumulative 
squeezing is the only one (or the most effective) squeezing 
strategy. We have studied the best localization that can 
be achieved with a given number of identical (5-kicks, by 
minimizing the localization factor Eq.(^l]) with respect 
to the delay times At„ between the kicks. For clarity, 
we present here only the results for zero initial tempera- 
ture of the atoms. Table 1 shows the best values of the 
localization factor found for up to five kicks, and com- 
pares them with the results of the accumulative squeez- 
ing strategy with the same number of kicks. While the 
maximal atomic localization that can be achieved with 
two pulses is almost the same for accumulative squeezing 
and for the optimal sequence of two pulses, the optimized 
results for three and more pulses are much better. 

For illustration, we choose the sequence of four optimal 
pulses to visualize the dynamics behind the localization 
process. Figure 6 (a) shows the spatial distribution of 
the ensemble of atoms at the time of arrival of the sec- 
ond pulse. Note that the second pulse is not applied 
at the time of the maximal localization. On the con- 
trary, the optimized procedure finds it favorable to wait 
after the focusing event, until the distribution becomes 
rather broad. The optimal four-pulse sequence we found 
requires applying the third and the forth pulses simulta- 
neously, thus producing an effective " double pulse" . The 
spatial distribution at the time of the combined third and 
fourth pulses can be seen in Figure 6 (b). The distribu- 
tion is only slightly more localized than in Fig. 6(a), and 
only the last pulse (with the double strength) squeezes 
the ensemble at the time of the maximal localization, 
thus bringing most of the atoms to the optical lattice 
minima [see Fig. 6 (c)]. 

IV. ORIENTATION AND ALIGNMENT OF A 3D 
ROTOR 

Under certain conditions, the process of molecular ori- 
entation (or alignment) by laser fields can be described 
by a strongly driven three-dimensional rigid rotor model. 
Although many features in the dynamics of 3D rotors are 
similar to those already discussed for the two-dimensional 
case, there are two principal differences that we want to 
emphasize. The first one may be traced in the evolution 
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of a classical ensemble of 3D-rotors being initially at zero 
temperature. The probability of finding a rotor (driven 
by a linearly polarized field) at a certain solid angle el- 
ement sinOdO is determined by the initial distribution 
function, fo{do) as follows: 



fi9) smdd9^ foieo)sm9od9o. 



(13) 



Therefore, the probability density of finding a rotor at 
the angle 9 at some latter time is 



m = E 



./o(gg)singg 
\d9/d9§\ sin 9 



(14) 



where the summation is done over all branches of the 
9o{9) function (see Section II). The probability density 
f{9) has a singularity if one of the factors in the denomi- 
nator is zero. As in the 2D case, the zeroes of the fraction 
\d9/d9Q\ give rise to the formation of the rainbow-like 
structures. But in addition, the geometrical factor sin0 
can get zero too, which leads to the additional singular- 
ities at = and 9 = 7:. This kind of singularities are 
responsible for the formation of the corona and glory ef- 
fects in the optical and quantum-mechanical scattering 

The second difference appears at finite initial temper- 
ature of the ensemble of 3D rotors. Because of the con- 
servation of the angular momentum projection onto the 
field polarization direction, an effective repulsive centrifu- 
gal force appears that prevents the rotors from reaching 
the exact 9 = Q and 9 = n orientations. As a result, two 
holes in the angular distribution of a driven 3D thermal 
ensemble should be always present &l 9 = and 9 = tt 

il- 

We present below the results of Monte-Carlo simula- 
tion of the dynamics of a classical ensemble of three- 
dimensional rotors (linear molecules) driven by a lin- 
early polarized time-dependent field. The corresponding 
Hamiltonian is 



H = 



1 




— fi£(t)r cos9, 



(15) 



where m is the reduced mass of the molecule, r is 
the (fixed) distance between the atoms, fi is permanent 
dipole moment, 9 and are Euler angles, and pg and 
Pcj, are related canonical momenta. For the driving field 
described by a 5— pulse, the equations of motion can be 
easily integrated (see [|5j). 

At zero initial temperature, pg^, = and = 0, and 
the angle (f)(t) is an additional constant of motion. In this 
case, the dynamics of the system is reduced to that of a 
two-dimensional rotor, beside the effect of the geomet- 
rical factor described above. Figure 7 demonstrates the 
time evolution of the angular distribution of an ensemble 
with finite initial temperature. Here, the probability of 
finding a rotor inside the solid angle element sin 9d9d(j) is 



plotted on a sphere. Red color corresponds to the high 
probability density while the blue color presents a low 
probability. Initial isotropic ensemble corresponds to an 
uniform solid angle distribution [see Fig. 7 (a)]. In Fig- 
ure 7 (b), the distribution is plotted at the "focal time" 
defined according to Eq.(^). In the two-dimensional case 
considered in the previous Sections, the distribution at 
the focal time is characterized by a sharp peak at 9 = 
and a broad background for larger values oi 9. In the 
three-dimensional case, the distribution has an additional 
sharp dip at 9 = 0, as discussed above. Figures 7 (c) and 
(d) show the angular distribution at t > tf. In Fig. 7 
(c) the formation of the " corona!^ around 9 = can be 
seen. In addition, we can observe a ring of relatively high 
probability moving from the north pole to the south pole 
of the sphere. This ring is a three-dimensional analog 
of the rainbow structure modified by the thermal effects. 
After the ring arrives at the south pole, a singular fea- 
ture appears around 9 = tt, with a hole in the center 
caused the repulsive centrifugal force [Fig. 7 (d)]. The 
formation of this sharp and robust structure in the an- 
gular distribution is analogous to the glory effect in the 
wave optics. 

All the described phenomena, namely, focusing, caus- 
tics creation, and accumulating squeezing are not specific 
for the simplest models considered above, but are rather 
common features that can be observed under general con- 
ditions of a strong excitation of the quantum rotor. In- 
deed, for strong enough driving field one may neglect the 
initial rotational energy stored in the rotor, and use the 
above quasiclassical ideas to describe its dynamics. The 
role of the initial state of the rotor (even thermally av- 
eraged) is reduced to the formation of a frozen classical- 
like initial angular distribution of the rotational ensem- 
ble. For example, for the dipole-type interaction, Eq. 
(^), the angular focusing may be analyzed by consider- 
ing classical dynamics at small angles {9 ^ 1): 



■e(r) 



0. 



(16) 



Here, r = th/I, and e = fi£(t)I /h^ , and the initial condi- 
tion is 9{tq) = 9q, d9{TQ)/dT = (where tq is any moment 
in the past before the beginning of the pulse). Focusing 
time Tf is defined by 9{Tf) = 0. Because of the linear- 
ity of Eq. (p^ , the position of the focusing times, and 
number of focusing events do not depend on 9o, but are 
determined only by the properties of the pulse (.(t). This 
boundary problem may be solved analytically in several 
special cases only (for (5-pulse considered above, or for 
a step-like e(T) dependence). In general, determination 
of focusing times requires numerical solution of Eq. (^|) . 
Based on the previous analysis, we expect that a new 
angular rainbow appears immediately after each focus- 
ing event. 

Figure 8 shows numerically calculated time-evolution 
of a three-dimensional quantum rotor (linear molecule) 
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coupled to the external laser field via the anisotropic po- 
larization interaction. The time averaged value of this 
interaction is (see, e. g. ]2l|J^) 

V{e,t) = -l/4£2(t)[(Q,|| -ajL)cos^{0) +ai_] (17) 

Here a|| and a±_ are the components of the polarizabil- 
ity, parallel and perpendicular to the molecular axis, and 
£{t) is the envelope of the laser pulse. In this case there 
arc two opposite equilibrium directions 6 ~ 0,tt, so the 
field aligns the molecules but not orients them. The ap- 
plied pulse has a Gaussian shape of a finite duration (i. 
e. it is not a (5-kick). As seen from Fig. 8, both focusing 
and caustics formation occur in the system, similar to the 
previous consideration, with the differences attributed to 
another symmetry of the problem. For instance, the an- 
gular rainbows are made of moving rings located sym- 
metrically with respect to the equatorial plane. Note, 
that the glory and corona are not seen because of the 
sin(0) factor incorporated into the distribution function 
shown in Fig. 8. 



V. CONCLUSIONS 

The predicted effects may be observed in a wide range 
of systems with strongly driven rotational degrees of free- 
dom. Possible examples range from heavy-ion collisions 
(when highly excited wave packets of nuclear rotational 
states are produced to molecules subject to strong 
laser pulses, and cold atoms trapped by standing light 
waves. The spectacular features described in this pa- 
per may be observed in the spatial distribution of an 
atomic ensemble driven by pulsed optical lattices. Re- 
cently, the accumulated squeezing scenario fl^ has been 
realized experimentally in this system [Q. Moreover, 
the related squeezing approaches may find application in 
atom lithography of ultra-high resolution [^,^ . In the 
case of molecules, the considered effects may reveal them- 
selves in the angular distribution of fragments produced 
by intense laser-field molecular interaction. The most di- 
rect evidence can be achieved in a two-pulse experiment, 
in which the first strong non-resonant pulse attempts to 
orient (align) the molecular ensemble, while the second 
short delayed pulse creates fragment ions. 
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Figure Captions: 

Table 1. The table shows the minimal value of the lo- 
calization factor that can be achieved with a fixed num- 
ber of kicks using accumulative squeezing scenario (Oacc) 
and optimized sequence of pulses (Oopt)- 

Fig. 1. Angular distribution of a quantum rotor ex- 
cited by a strong ^-kick (P = 85). The graphs corre- 
spond to (a) T — O.bTf, (b) T — Tf, (c) r = 2t/, (d) 

T = Tf+Trevf^, (c) T = Tf+Trev/i, (f) T = Tf+Trev/'i, 

(g) T = 1.8t/ + Trev/2, (h) T = I.Stj + Trev/S, and (i) 
T — I.Stj -|-Tret,/4, respectively. 

Fig. 2. Classical map representing the final angle 6 as 
a function of initial angle for (a) r = 0.5ry, (b)r = t/, 

(c) T = 3t/, and (d) r = IOtj. 

Fig. 3. Accumulative angular squeezing. Graphs are 
shown in double logarithmic scale. 

Fig. 4. Spatial distribution of a classical ensemble of 
atoms after a single ^-kick. In Figs, (a) and (b), the 
initial temperature, T of the ensemble is zero. In (c) and 

(d) , keT = il/9)kf[J^^V{t)dt]^/m, where V{t) is the 
depth of the potential produced by the standing wave, 
and fcs is Boltzmann constant. Figures (a) and (c) show 
the spatial distribution at focal time r = t/, while in 
figure (b) and (c) r = 2.5r/. 

Fig. 5. Accumulative squeezing of atoms in a 
pulsed optical lattice at finite temperature (classical de- 
scription). The minimal localization factor as a func- 
tion of the kick number is shown in double logarith- 
mic scale. The solid line corresponds to zero ini- 
tial temperature. The dashed and dotted lines corre- 
spond to keT = {l/9)kf[J^^V{t)dt]'^/m and keT = 

(4/9)fc?[/^ Vit)dt\^/m, respectively. 

Fig. 6. Spatial distribution for the optimized sequence 
of four 5-pulses. The upper row shows the spatial dis- 
tribution averaged over 100 atomic ensembles (each con- 
taining 5000 atoms), the lower figures show the distribu- 
tion of atoms in one of the ensembles. In (a), the dis- 
tribution is plotted at the time of the second pulse, that 
is delayed by r = 3.02r/ with respect to the first pulse, 
in (b) - at the time of the (combined) third and fourth 
pulses delayed by r = 1.35r/ after the second pulse. Fig- 
ure (c) shows the distribution of atoms at the time of the 
maximal squeezing. 



Fig. 7. Time evolution of the angular distribution 
of a (classical) ensemble of (5-kicked 3D rotors at finite 
temperature. The probability of finding a rotor inside 
the solid angle element sv[\6d6d(j) is plotted on a sphere. 
Picture (a) shows the distribution at r = (time of the 
kick). Figures (b), (c), and (d) correspond to r = t/, r = 
3.3Ty and r ~ 5Tf, respectively. The initial temperature 
corresponds to ksT = {1 / 100) n^r^[ J ^^£{t)dt]^ /P. 

Fig. 8. Contour plot for the time-dependent 
angular distribution function 27rsin(6')|^'(0, r)|^ of a 
three-dimensional rotor (molecule) subject to a strong 
" polarization- type" interaction, Eq.(|l7|) with a Gaus- 
sian pulsed laser field: £'^{t){a\\ — a±)I/Ah^ = 3 x 
10^ exp[— (r/0.01)^]. Here t — tti/I. The molecule re- 
sides initially in the isotropic ground angular state (J = 
0, m = 0). Angular focusing and rainbows emerging from 
each of the focal regions can be seen. The angular rain- 
bows are made of moving rings located symmetrically 
with respect to the equatorial plane. 
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No of kicks 


Oacc 


Oopt 


2 


0.33 


0.31 


3 


0.26 


0.20 


4 


0.21 


0.11 


5 


0.18 


0.07 



TABLE I. 
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